Effects of Matter Density Fluctuation in Long 
Baseline Neutrino Oscillation Experiments 

oo ; Masafumi Koike a *and Joe Sato &t 

ON ■ ° Institute for Cosmic Ray Research, University of Tokyo, 



OO 



X 



Midori-cho, Tanashi, Tokyo 188-8502, Japan 
Department of Physics, University of Toky~ 
Bunkyo-ku, Hongo, Tokyo 133-0033, Japan 



b Department of Physics, University of Tokyo, 



> 

y—i . Abstract 
<N 

£f) ' The effects of matter density fluctuation in long baseline neutrino oscillation 

experiments are studied. Effects of short wavelength fluctuations are in general 
irrelevant. Effects of long wavelength fluctuations must be checked on a case-by-case 
basis. As an example we checked the fluctuation effects and showed its irrelevance 
in a case of K2K experiments. 
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Neutrino oscillation search is one of the most promising way to go beyond the 
Standard Model. The existence of neutrino oscillation is suggested by the lack of 
solar neutrino flux [jy, ^, £3, |j and the anomalies of atmospheric neutrinos [||, 
[/J []. Further confirmation of its existence is now close at hand owing to the long 



baseline neutrino oscillation experiments, such as K2K experiments [1C] and Minos 
experiments pl| . The possible effects of CP and T violation in such experiments 
have been discussed in the literatures |[^, 13, |l4|, [H], [l^, 17 1. There the matter 



density is approximated to be constant^. We now study effects of the matter density 
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2 We have checked numerically in Ref. |17| that the constant density approximation works well for 
some parameter pace of neutrino's in K2K experiments. However it does not mean the approximation 
works well for all the parameter space. Another example of numerical calculations is presented in Ref. 

in. 
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fluctuation and clarify when the constant density approximation is valid. We also 
verify that this approximation works well for the K2K experiment. 

Let v = (i/ e , u T ) and v 1 = {y\, V2, ^3) be the eigenstates of flavor and mass in 
vacuum, respectively. They are related through the mixing matrix U as 



v = Uv'. 



(1) 



Defining Srrifj = mf — mh with the mass eigenvalues rrii(i = 1,2,3), the evolution 
equation of flavor eigenstates in matter is given by 

. dv(x) 



dx 



H{x)v{x), 
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H{x) = 
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5m 2 



21 



E7+ + 



Here 



a(x) = 2V2G F n e (x)E = 7.56 x 10~ 5 eV 2 



a(x) 



p(x) 



E 



(2) 
(3) 

(4) 



g cm 3 GeV ' 

n e (x) is the electron density, p(x) is the matter density and E is the neutrino energy. 
The evolution equation @ is solved as 



v{V) = Texp ^-i dxH(x)j z/(0) = S(L)v(Q), 



giving the oscillation probability from v a to vp(a,P = e,fi,r) as 

P{y a -> vp; L, E) = \S(L)p a \ 2 . 



(5) 



(6) 



We assume Sml x ~ (10~ 2 ~ l(T 3 )eV 2 and bm\ x ~ (10" 5 ~ 10~ 4 )eV 2 al- 
lowing for atmospheric neutrino anomaly and solar neutrino deficit. Along with 
a ~ 10~ 4 eV 2 (see eq.@), we see 

Sm^i, a <C Sm 3 i. 

We hence separate as H(x) = Hq + H\{x), where 

, / 



(7) 



H = — U 
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and 



H x {x) = 



2E 



U 



( 



V 



5m 2 
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(8) 



(9) 



and treat H\[x) as a perturbation. Taking up to the lowest order in H±(x), we 
obtain 



S(L) 



~ e~ iHoL + e~ iHoL (-i) 



5 (L)+5i(L). 



dx 



e iHoX Hi(x)e- iHoX 



(10) 
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The fluctuation of matter density p(x), or equivalently that of a(x), affects S\(L) 
alone. We separate that effect from Si(L) in the following. 
Letting 

1 f L 

a = — / dx a(x) 
L Jo 



and 

Sa(x) = a(x) — a, 
we separate the contribution of 5a(x) from H\{x) by defining 

/ 



Hi(x) 
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and 



Defining accordingly 



8H 1 {x) = 



2E 



Si 



(L) = -i f 
Jo 



dx 







jf/t+|^ a o 








'•.)■ 








- x) ^i(x)e- iHox " 





and 



SSi(L) = -i / dx 
J o 

Si is calculated to be Q| 

Sl{x)f3a = -iU{3 2 U* 2 - 



e -iH (L-x) SHi ^ e -iH 0a 



8m%xX 



2E 



-i(5p e - Ufs 3 U* 3 )(8 ae - U e3 U* 3 ) 



ax 
2E 



exp — i 
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+U f33 U: 3 (5 ae -U e3 U* a3 ) 
+(^e - Up 3 U* c3 )U c3 W a3 

To evaluate 8S\(x)p a we expand <5a(x) as 

oo 

<5a(x) = a ne~ innx/L 



5m 31 x^ 


- 1 
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(11) 

(12) 
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(15) 
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and carry out integration for each Fourier component. Note that ao = from the 
definition (12). We then find 



SSi(x)p c 



8aeU/33U* 3 



a 3 



cxp 



cxp 



U{3sU* 3 \ U e3 \ 



exp 



. 5m 31 x 
1 2E 
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(19) 



We now see the order of magnitude of S\ and <5Si(x) ((0) and (|l~9|) ). Here we 
assume that the products of t/'s and dm^L j2E are all 0(1) so that we can observe 
the neutrino oscillation. In this case we can see from (|l^) that 



S 1 = 



5m? 



21 



Smly 



or 



(20) 



On the other hand all the three terms of SS\(x) contains factors (a n /5m 3 i) and also 
{1 - 2mr/(5ml 1 /2E)}- 1 ~ 1/n. Hence 



n a 



(21) 



We can see from (pl| ) that the long wavelength fluctuation (i.e. small n) of the 
matter density is important; short wavelength fluctuations (i.e. large n) of matter 
density is in general irrelevant due to the factor 1/n. This means we do not need 
to survey detailed profile of the matter density distribution on the baseline. On 
the other hand we must check on a case-by-case basis whether the long wavelength 
fluctuation is relevant or not. This check can be done by considering the magnitudes 
of a n /a. 

Let us carry out a check for the K2K experiments as an example. Figure |l| 



shows the density profile between KEK and Kamioka [19]. The first several Fourier 
coefficients divided by the mean density (= a n /a) for this profile is given in Fig. 
^. We can see that a n /a is small, mostly much less than 0.1. This justifies the 
constant matter density approximation for K2K experiments. 

We have shown that short wavelength fluctuation of matter density is irrelevant 
and only first several Fourier coefficients of fluctuation may be important for long 
baseline neutrino oscillation experiments. Since we can expect that such fluctuations 
are much smaller than the mean density in crust we can approximate the matter 
density to be constant. This is indeed the case for K2K experiments. We could check 
the approximation works well independent of the parameters relevant to neutrino 
oscillation. 
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Figure 1: The density profile between KEK and Kamioka. 











• 
















• 

• 








• 1 


------- 







10 



15 



20 



Mode Index n 



Figure 2: The values of a n /a for first several n's. 
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